Superconductivity in repulsive electron systems having three-dimensional 

disconnected Fermi surfaces 



m 
o 
o 



o 
o 

Oh' 



o 
o 



> 

(N 

m 
O 
m 
o 

-I— > 
c3 



i 

C 

O 

o 



Seiichiro Onari 1 , Ryotaro Arita 1 , Kazuhiko Kuroki 2 , and Hideo Aoki 1 
1 Department of Physics, University of Tokyo, Hongo, Tokyo 113-0033, Japan 
2 Department of Applied Physics and Chemistry, University of Electro-Communications, Chofu, Tokyo 182-8585, Japan 

(February 2, 2008) 

The idea of raising T c in the spin-fluctuation mediated superconductivity on disconnected Fermi 
surfaces with the gap function changing sign across but not within the Fermi pockets, proposed by 
Kuroki and Arita for two dimensions (2D), is here extended to three-dimensional (3D) systems. Two 
typical cases of 3D disconnected Fermi surfaces (stacked bond-alternating lattice and stacked ladder 
layers) are considered. By solving Eliashberg's equation for Green's function obtained with the 
fluctuation exchange approximation (FLEX) for the repulsive Hubbard model on these structures, 
we have shown that T c can indeed reach O(0. Olt), which is almost an order of magnitude higher 
than in ordinary 3D cases and similar to those for the best case found in 2D. The key factor found 
here for the favorable condition for the superconductivity on disconnected Fermi surfaces is that 

the system should be quasi-low dimensional, and the peak in the spin susceptibility should be 
appropriately "blurred" . 

PACS numbers: 74.20.Mn 



I. INTRODUCTION 



The monumental discovery of high-T c superconduc- 
tivity in strongly correlated electrons in cuprates 1 has 
kicked off intensive studies on superconductivity in repul- 
sively interacting electron systems. These studies have 
brought to light the spin-fluctuation mediated pairing 
mechanism, which gives rise to anisotropic Cooper pair- 
ing instabilities. So identifying the optimum condition 
for pairing in the repulsively interacting electron systems 
is one of the most fascinating goals of theoretical studies, 
since it should give a guiding principle in searching for 
higher- T c materials. 

A fundamental question in this context is a repeatedly 
asked one: "why is T c in the electron-mechanism super- 
conductivity so low?" Although this may at first sound 
strange, this is true in that the estimated T c for ordi- 
nary lattices is at best O(0. Olt) 2 which is two orders of 
magnitude smaller than the starting electronic (hopping) 
energy, t. There are in fact good reasons why T c is so low. 
One is that the spin-fluctuation mediated interaction is 
much weaker than the original electron-electron interac- 
tion. Another important reason, also quite inherent in 
the electron mechanism, is the following. If one wants to 
have superconductivity from repulsion, pairing has to be 
anisotropic, since the nodes in the gap function, which 
convert the repulsion to an attraction in the gap equa- 
tion, realize the superconductivity. Namely, the pair- 
scattering processes (the matrix elements of the inter- 
action that connect two electrons around the Fermi sur- 
face to another positions) contribute positively to the gap 
equation when there is a node across the pair-scattering 
momentum transfer. Unfortunately, the nodes greatly re- 
duce T c : While the gap function has opposite signs across 
the main pair scattering (e.g., k — (0,7r) <-» (ir,0) for 



the antifcrromagnetic spin- fluctuation exchange), there 
are other pair scattering processes across the regions of 
the Fermi surface having the same sign, and these work 
against the pairing and lowers T c . Since the nodes are 
imperative in anisotropic-pairing superconductivity, this 
problem is unavoidable as far as the ordinary lattices are 
concerned. 

So a novel avenue to explore is: can we improve the 
situation by going over to multiband systems. In this 
context Kuroki and Arita 3 have proposed that the sys- 
tems that have disconnected Fermi surfaces are a promis- 
ing candidate for higher- T c superconductivity. In such 
systems the gap function can change sign across the 
Fermi surfaces but not within each Fermi surface, so 
that all of the main pair-scattering processes contribute 
positively to superconductivity. The estimated T c for 
two-dimensional (2D) repulsive Hubbard model on such 
lattices is indeed almost an order of magnitude higher, 
T c ~ O.li. The disconnected Fermi surface has been sub- 
sequently explored in other models and systems 4-6 . 

As for the dimensionality of the system, on the other- 
hand, three of the present authors have systematically 
studied 2 the Hubbard model on 2D and 3D single- 
band lattices with the fluctuation-exchange approxima- 
tion (FLEX), and have clarified that 2D systems are 
generally more favorable than 3D systems as far as 
the spin-fluctuation-mediated superconductivity in ordi- 
nary lattices (square, triangular, fee, bec, etc) are con- 
cerned. The "best case" among the ordinary lattices is 
identified to be the d-wave pairing in the square lat- 
tice, and the above-mentioned T c ~ O(0. Olt) in fact 
refers to this case. Physical reason why 2D is more 
favorable is that the fraction of the phase volume over 
which the fluctuation-mediated interaction exerts signifi- 
cant attractions is much smaller in 3D. 2 T c estimated by 
solving Eliashberg's equation with FLEX is indeed only 
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O(O.OOli) even for the best 3D case (as far as simple 
lattices are concerned). Similar results have also been 
obtained by Monthoux et al. 7 independently by means 
of a phcnomcnological calculation. Recent experimen- 
tal finding 8 that a series of heavy-fermion compounds 
Ce(Rh,Ir,Co)In5 has higher T c when the structure is more 
two-dimensional (with larger c/a) is consistent with this 
prediction (although orbital degeneracies in the heavy- 
fermion compound may exert some effects). 

Now, if one puts these two observations (i.e., one on 
the disconnected Fermi surface and another on the di- 
mensionality), a natural question arises: can we con- 
ceive 3D lattices having disconnected Fermi surfaces that 
have high T c 's. In other words, can the disconnected 
Fermi surface overcome the disadvantage of 3D to possi- 
bly realize T c ~ O(0. It)), which should be a fundamental 
question for the spin-fluctuation mediated pairing. More 
specifically, what we have in mind is the inter-band nest- 
ing in the 3D disconnected Fermi surface as schematically 
depicted in Fig.l. 




kx 

FIG. 1. Schematic inter-band nesting (arrows) in discon- 
nected Fermi surfaces in 3D. 



There, a nesting vector exists across the two bands, 
and this is envisaged to give the attractive pair-scattering 
interaction with a nodal plane running in between the 
two bands. Note that a pair is formed within each band, 
although the nesting is inter-band. In such a case the 
gap function should be nodeless within each band, while 
the gap has opposite signs between the two bands. In 
this sense the mechanism is an extension of the time- 
honored Suhl-Kondo mechanism, although there are im- 
portant differences since here we consider bands originat- 
ing from multiple orbitals of the same type connected by 
single electron hopping in order to have appropriate nest- 
ing between the disconnected Fermi surfaces, while Suhl- 
Kondo considers basically the bands originating from dif- 
ferent (e.g., s and d) orbitals that interact only via multi 
electron processes. 

For the survey in lattice structures, there are two pos- 
sible approaches: (i) to find 3D lattice structures whose 
Fermi surfaces are disconnected, and consider the 3D- 
version of the mechanism originally proposed for 2D 
lattices, 3 (ii) stack the 2D lattices already proposed 3-6 
for which T c ~ 0(0. It) has been estimated, and intro- 
duce the interlayer hopping, t z , to see how the T c be- 



haves as t z is increased. As for the former approach, we 
have recently studied superconductivity in 3D stacked 
honeycomb lattice, 9 where the Fermi surface consists of 
two networks of tubes. However, the networks of tubes 
have turned out to have such a strong structure that 
the intra-band as well as inter-band pair-scattering pro- 
cesses contribute to the gap equation, and the resulting 
gap function has nodes (with d-wave symmetry) within 
each band. Accordingly, the estimated T c ~ O(0.001i) is 
low, and the stacked honeycomb lattice does not realize 
the 3D-version of inter-band nesting conceived here. 

This has motivated us in the present study to look for 
models that have more compact and structure-less Fermi 
surfaces. We have found that a stacked bond-alternating 
lattice (Fig. 2) has a compact and disconnected Fermi 
surface (i.e., a pair of ball- like Fermi pockets). We then 
show for the multiband Hubbard model on this lattice 
that the inter-band pair-scattering alone is dominant. 
The estimated T c is 0(0. Olt), which is the same order 
of that for the square lattice, and remarkably high for a 
3D system. 

As for the second approach of stacking the 2D lat- 
tices, we have studied the ^-dependence of the pairing 
instability for stacked ladder lattice (Fig. 2). We have 
found that the T c ~ 0(0. If), found for the single layer, 
is unexpectedly robust against, or even increases with, 
t z when t z is not too strong. We have then identified 
the underlying physics for the favorable condition for the 
superconductivity on disconnected Fermi surfaces — the 
system should be quasi-low dimensional, and the peak in 
the spin susceptibility should be appropriately "blurred" . 

So we shall conclude that the idea of raising T c in su- 
perconductivity from repulsive electron-electron interac- 
tion on disconnected Fermi surfaces can work in 3D lat- 
tices if we consider appropriate lattice structures. 



II. FORMULATION 



Here we take the repulsive Hubbard model, a simplest 
possible model for repulsively interacting electron sys- 
tems. So the interaction is the ordinary on-site repul- 
sion, £7, but the lattice has two atoms per unit cell here 
accompanied by two bands. The Hamiltonian is 

nn 

« = EE E 

i,j cr=T,i a,/3=A,B 

+ U J2 E <T<i> (!) 

i a=A,B 

where a, (3 label A,B sublattices with standard notations 
otherwise. For the transfer energies t, t Xl , t X2 and t y are 
intra-layer hoppings, while t z the inter-layer (see Fig. 2). 
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FIG. 2. A stacked bond-alternating lattice (left panel) and 
a stacked connected-ladder lattice (right) considered in the 
present work. Solid lines (dashed lines) indicate stronger 
(weaker) hopping, and A, B label sublattices. 



We employ the FLEX method developed by Bick- 
ers et al. 10 ~ 13 In the two-band version of FLEX, 14 ' 15 
Green's function G, self-energy S, susceptibility x, and 
the gap function <p ah become 2x2 matrices, such as 
G a p(k), where a, (3 = A,B sublattice, k = (k,iiv n ) with 
w n = (2n — l)nT being the Matsubara frequency for 
fermions. 

The self-energy is given in the FLEX as 



-V pp (k-k')Gp a (-k'] 



(2) 



where the particle-hole (V ph ) and particle-particle (V pp ) 
fluctuation-exchange interactions are given as 



X ph (q) 



1 - U X ^{q) 



1 , 



a/3 



X ph (q) 



1 + U X ph (q) 



-u 2 x%(q) 



a/3 

(3) 



with the energy of a free electron. 

The gap function, <fr, and T c are obtained from Eliash- 
berg's equation, 

= EE 



N 



k> a', 13' 



V${k - k')G aa/ (k')G^(-k')<p a , ,(k'), (9) 

where the (spin-singlet) pairing interaction V^(q) is 
given as 



X ph (q) 



1 - U X ph (q) 



1 



a/3 



x ph (q) 



l + u x ph (q) 



a/3 

(10) 



T c is then determined as the temperature at which the 
largest eigenvalue, A, of Eliashberg's equation becomes 
unity. 

The spin susceptibility in RPA is given by, 



x ph (fc,Q) 

l-U X ph (k,0) 



(11) 



a/3 



may be expressed as diagonalized components, 
Xaa + Xbb 



X± 



± 



XAA - XBB 



+ lx 



AB | 



(12) 



For the two-band FLEX calculation in 3D we take N = 
32 3 fe-point meshes, and the Matsubara frequencies taken 
as -(2N C - l)nT < iv n < (2N C - 1)ttT with iV c = 4096, 
which gives converged results. 



V P ${q) = U 1 



X pp (q) 



1 + U X pp {q) 



X PP {q) 



a(3 



with 



xtpiq) 



-J2G al3 (k + q)G 0a (k), 



T 



^(9) = ^E G ^(fc + <z)G Q/3 (-fc)- 



(4) 

(5) 
(6) 



Here we denote q = (q,iei) with e; = 27r/T being the 
Matsubara frequency for bosons, and N the number of 
fc-points on a mesh. 

Green's function is related to S via Dyson's equation, 



[G(k)- 1 ] a/3 = [G°(k)- 1 ] af3 -i: a0 (k), 

where G° is the bare Green's function, 

1 



G%(k) = 



iuj n + /j, - e 



(7) 



(8) 



III. RESULT 

A. Stacked bond-alternating lattice 

We have searched for the case where the inter-band 
nesting is strong with a weak intra-band nesting along 
the following line. We start from the stacked honey- 
comb lattice, where the tubular Fermi surface connected 
along x — y directions has too strong an intra-band nest- 
ing as mentioned. 9 One way to make the Fermi surface 
more compact and structureless is to introduce alternat- 
ing bonding structures in the x — y plane. For example, if 
the in-plane bonds along x alternate between weak ones 
{t x \) and strong ones {t X 2), where the array of strong 
and weak ones is staggered as we go along y direction. 
We can then connect the layers with the third-direction 
hopping t z , and we call this the stacked bond-alternating 
lattice (Fig. 2). We have varied t Xl and t z to examine the 
pairing instability over the region t Xl ,t z < t y (= 1) < t X2 
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with a fixed t X2 = 2 for the band filling n = 1.02. This re- 
gion has turned out to correspond to the case of compact 
Fermi pockets. The repulsion is taken to be U — 4, which 
is a typical value for strongly correlated systems. Here- 
after we take t y as the unit of energy. The Fermi surface 
and the spin susceptibility for a typical case (t Xl = 0.4, 
t z = 0.8) is shown in Fig. 3. We can see that the Fermi 
surface comprises a pair of compact, ball-like pockets. 
There exists, as intended, a dominant nesting vector in 
the stacking direction (k = (0, 0, tt)) across the two pock- 
ets, as confirmed from a strong peak in X+(fe, 0) around 
k = (0,0, tt) in Fig. 3. 



though disconnected Fermi surfaces are favorable for the 
inter-band pairing, the Fermi pockets have to be compact 
so that the intra-band pair-scattering is suppressed and 
the gap function for each band becomes nodeless. 



7=0.02 U=4 72=1.02 ta=0.2 ta=2.0 




FIG. 3. (Color) Fermi surface (left panel; blue: bonding 
band, red: antibonding band) along with x+ at k z = n (right) 
for the stacked bond-alternating lattice with U = 4, n = 1.02, 
t xi = 0.4, t xa = 2.0, t z = 0.8, T = 0.02. 



Figure 4 shows the result for the dependence of the 
largest eigenvalue, A, of Eliashberg's equation on the 
inter-layer hopping, t X) for T = 0.02 and t xi = 0.2. We 
can see that A increases with t z . This indeed occurs as the 
two Fermi pockets (displayed in red and blue) begin to 
be formed and the inter-band nesting across the pockets, 
that should favor the pairing, develops. For this choice 
of parameters, however, the superconductivity gives way 
to antiferromagnetism (whose boundary is identified here 
to the situation when 1 — Ux ph = 0) around the value 
of t z at which A reaches unity. Nevertheless we are on 
the right track in that the gap function is, as intended, 
nodeless (gapful) in each band while changes sign across 
the two bands. 

So we should next tune the shape of the two Fermi 
pockets. This can be achieved by varying the weak hop- 
ping, t xi , in the bond-alternating system. As seen in Fig. 
5, the Fermi surface for each band remains to be two balls 
throughout, but they become more compact as t Xx is in- 
creased. If we turn to A in the figure for a fixed t z = 0.8 
with T = 0.02, A indeed reaches unity well before the 
antiferromagnetism sets in. In the case of t x x = (with 
a low T c ), the gap function has a nodal plane as shown 
by green sheet in Fig. 5, which should be due to the 
main (0, 0, tt) pair scattering in large antibonding band 
Fermi surface favor node to change the sign across the 
scattering. 

So we conclude for the bond-alternating model that, al- 




f--=0-2 fc=o.6 f 2=0.8 

FIG. 4. (Color) The largest eigenvalue, A (top left), of 
Eliashberg's equation and the peak value of spin susceptibility 
X+ (top right, with the gray region indicating the antiferro- 
magnetic phase) as a function of the inter-layer hopping t z , 
along with the Fermi surface (bottom) for selected values of 



t z for U = 4, n = 1.02, U 



0.2, t x 



2.0, T = 0.02 for the 



stacked bond-alternating lattice. 
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r=0.02 U=A n=1.02 &2=2.0 fe=0.g 




o fcv 

1x1=0.2 

FIG. 5. (Color) A plot similar to Fig. 4 (top two 
rows) when we vary the weak hopping, t xi , in the stacked 
bond-alternating lattice with fixed t X2 = 2.0, t z = 0.8 for 
U = 4, n = 1.02, and T = 0.02. We depict the gap function 
with its sign indicated on the Fermi surface for the bonding 
(second from bottom row) and antibonding (bottom) bands, 
respectively, where nodal planes are shown by green sheets. 



We have then scanned the t xi — t z parameter plane to 
optimize the values of t Xl and t z , and found that A takes 
its maximum at t Xl — 0.4, t z = 0.8. In Fig. 6, temper- 
ature dependence of A for the optimized case is plotted. 
We can see that T c ~ 0.02^, which is remarkably high 
(i.e., ~ 2 x 10 _3 x the bandwidth with the bandwidth 
~ 12t y ) for a 3D system (where usually T c ~ 10 _4 x 
bandwidth at most as estimated for the simple cubic 
lattice 2 ). 



U=4 n=1.02 fci=0.4 ta=2,0 fe=0-8 




FIG. 6. A as a function of T for the optimized t xi — 0.4 and 
t z = 0.8 for the stacked bond-alternating lattice with [7 = 4, 
n = 1.02, tx, = 2.0, and T = 0.02. 



B. Stacked connected-ladder lattice 

Let us move on to the stacked connected-ladder lat- 
tice (right panel in Fig. 2). In this subsection, we fix 
n = 1.05, U = 4.0, t X2 = 1.6, tj, = l(unit). This value of 
t X2 turns out to give the highest T c for U — 4. In 2D 
Kuroki and Arita 3 have shown that a larger U — 8 gives 
a higher T c , but in the present case of 3D the tendency 
toward AF is stronger, so that we take a weaker U = 4 
as in the previous section. 

Figure 7 shows A and \+ as a function of the inter- 
layer hopping t z for a fixed t x \ =0.1 with T = 0.06. 
We can see that A is as large as, or even greater than, 
in 2D (t z = 0) up to t z < 0.3. For smaller t z < 0.2 
the Fermi surface consists of only one band (two pan- 
els at bottom left in Fig. 7). In this case the inter-band 
nesting is weak, as seen from a small x+ (top right in 
Fig. 7). When t z > 0.3 ^> t x \, the system becomes quasi- 
2D along the yz plane, and a nesting in the direction of 
(0, 7T, 7t) evolves (as evident from the shape of the Fermi 
cylinders that change their direction from || z to || x). 
While x+ monotonically increases with t z throughout, 
superconductivity is most favored (A ~ 1 in between 
(t z ~ 0.2). 

The resultant T c ~ 0.06 in the 3D system is T c ~ 
8 x 10 _3 x bandwidth (with the bandwidth ~ 8t y ), which 
is even higher than in the above bond-alternating lattice, 
and almost an order of magnitude higher than the ordi- 
nary 3D case of ~ 10~ 4 x bandwidth. This shows that 
the superconductivity in the disconnected Fermi surface 
proposed in Ref. 3 for 2D is robust against introduction 
of the third-direction hopping. 
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C. Condition favoring the pairing on disconnected 
Fermi surfaces 

The peak in T c at a nonzero t z (rather than at t z = 0) 
seen above is surprising. This is particularly puzzling 
since the spin susceptibility increases monotonically with 
t z , which implies that a good nesting per se does not 
necessarily favor the pairing. We then want to pin-point 
the physical reason why superconductivity is degraded 
although x+ increases in the region of t z > 0.3. 

To elaborate the question let us look at another, similar 
phenomenon when we increase the in-plane, inter-ladder 
hopping, t xi , with a fixed t z = 0.2 in Fig. 8. There, 
we can see that, although the nesting, in the direction of 
(0, 7r, 7r), evolves with t Xl as shown from both the shape 
of the Fermi surface and the value of x+ i n Fig. 8, A de- 
creases as x+ increases, which provides another example 
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of the degraded superconductivity in the evolution of the 
nesting. 

We can conceive several reasons as follows. As we have 
mentioned, the increase of t z or t x i raises the dimension- 
ality of the anisotropic system to 3D. This should de- 
grade superconductivity in two ways. First one is just 
the mechanism revealed by Arita et at in ref. 2 , where 
higher the dimension smaller the phase volume fraction 
for the pairing interaction. Secondly, when Ep is close 
to band edges, the density of states there tends to be 
smaller for higher dimensions. To see that this is indeed 
the case, we have examined in Figs. 10, 11 the regions in 
fe-space in which 

pi - 8pL < e° + RcE < fi + 8fi, 

which we will refer to as the "thickness" of the Fermi sur- 
face. Bulkier the thickness larger the density of states, 
but we specifically take the energy window Sfj, as the en- 
ergy interval over which Imx + (fc max , ui) (with fc max being 
the momentum that gives the peak in x+)i displayed in 
Fig. 9, is significant. In other words, S/j, is the energy scale 
of the spin fluctuations that mediate the pairing. As seen 
in Figs. 10,11, respectively, increased t z or t x \ makes the 
Fermi surface "thinner" although the Fermi surface itself 
is large. The situation is shown schematically in Fig. 12 

Another factor governing superconductivity can also 
be identified from Figs.10,11. Namely, since we are here 
dealing with small (or vanishing) Fermi surfaces, the spin 
susceptibility is governed not by the shape of the Fermi 
surface itself, but by the thickened Fermi surface. The 
thin Fermi surfaces for smaller t z or t x \ then make the 
nesting sharper and the peak in x+ becomes narrower. 
This also acts to degrade superconductivity, since the 
spin susceptibility has to be appropriately "blurred" in 
order to make the spin fluctuation usable over a large 
phase space for pair scattering. 



FIG. 7. (Color) A (top left) and x+ (top right) as a function 
of the inter-layer hopping t z along with the Fermi surface for 
the stacked connected-ladder lattice with selected values of t z 
(as indicated at the bottom) for the stacked-ladder lattice for 
U = 4, n = 1.05, t xi = 0.1, t X2 = 1.6, T = 0.06. 



7"=0.1 U=4 n=\.Q5 te=1.6 tz=Q.2 
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tx\=0.\ m=0A hi=0.8 

FIG. 8. A plot similar to Fig. 7 for the stacked con- 
nected-ladder lattice when the weak hopping, t x i, is varied 
as indicated at the bottom for a fixed t z = 0.2 with U = 4, 



f/=4 H=1.02 ta=0.1 te=1.6 fe=0.2 



T=0.06 U-4 tl=1.05 txi =0.1 tx2=1.6 
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FIG. 9. Imx+(fc max , u) (normalized by its maximum 
value), where fe max = (0,7r,7r), as a function of ui for 
T = 0.06,17 = 4,w = 1.05, t Xl = 0.1,^ 2 = 1.6, t z = 0.2 
for the stacked connected-ladder lattice. 
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fc=0.0 tz=0.2 tz=0.4 

FIG. 10. (Color) Top: The regions in the fe-space bounded 
by fi—Sfj, < e°+RcE < fi+Sfj, with Sfi — 0.5(purple), l(green), 
along with the Fermi surface (blue: the bonding band which 
is only visible in the right row, red: the anti-bonding band) on 
which the region of e°+ReE = fi is displayed for the stacked 
connected-ladder lattice with selected values of the inter-layer 
hopping t z as in Fig.7 as indicated at the bottom. Bottom: 
fe-space contours (orange: \+ — 3, yellow: x+ — 1) of the 
spin susceptibility. 
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FIG. 11. (Color) A plot similar to Fig. 10 for the stacked 
connected-ladder lattice when the weak hopping, t xi , is varied 
as in Fig. 8 as indicated at the bottom. 
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FIG. 12. The schematic "thickness" of Fermi surface at 
k x — ir. The light gray and dark gray regions represent 
0.5 < |e° — fl + ReE| < 1 and < |e° - fi + ReE| < 0.5, re- 
spectively, and Fermi surfaces are shown by black lines. The 
left panel has larger gray (active) region than the right. 



IV. CONCLUSION 

We have studied how the pairing with higher T c in 
disconnected Fermi surfaces can be realized in three- 
dimensional systems. Our first proposal of the stacked 
bond-alternating lattice has a pair of ball-like Fermi 
pockets with a dominant inter-band nesting vector and 
antifcrromagnetic spin fluctuations, which is a genuinely 
3D model in that the interlayer hopping t z promotes 
the inter-band pairing. T c estimated from Eliashberg's 
equation shows that T c ~ 10~ 3 times the bandwidth 
(~ O(0. OH)), which is as high as that for the single-band 
Hubbard model on the square lattice, and remarkably 
high for 3D system. 

Our second model of the stacked layers each having a 
disconnected Fermi-surface in 2D originally proposed by 
Kuroki and Arita 3 is shown to have a high T c that is ro- 
bust against the introduction of the inter-layer hopping, 
although the hopping degrades the 2D nesting. The es- 
timated T c ~ several times 10 -3 x bandwidth (= several 
times O(0. Oli)) is even higher than the first model. 

From these we have then identified the right conditions 
for higher T c in the disconnected Fermi surface: the sys- 
tem should be quasi-low dimensional, and the peak in the 
spin susceptibility should be not too sharp but "blurred" . 
So the message here is that 3D materials with consider- 
ably high T c can be expected if we consider appropriate 
lattice structures. 
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